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Abstract 

Purpose Carbon footprint of field crops can be lowered 
through improved cropping practices. The objective of this 
study was to determine the carbon footprint of spring barley 
(Hordeum vulgare L.) in relation to various preceding oilseed 
crops that were fertilized at various rates of inorganic N the 
previous years. System boundary was from cradle-to-farm gate. 
Materials and methods Canola-quality mustard ( Brassica 
juncea L.), canola ( Brassica napus L.), sunflower ( Helianthus 
annuus L.), and flax ( Linum usitatissimum L.) were grown 
under the N fertilizer rates of 10, 30, 70, 90, 110, 150, and 
200 kg Nha -1 the previous year, and spring barley was grown 
on the field of standing oilseed stubble the following year. The 
study was conducted at six environmental sites; they were at 
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Indian Head in 2005, 2006 and 2007, and at Swift Current in 
2004, 2005 and 2006, Saskatchewan, Canada. 

Results and discussion On average, barley grown at humid 
Indian Head emitted greenhouse gases (GHGs) of 
1,003 kg C0 2 eq ha -1 , or 53% greater than that at the drier 
Swift Current site. Production and delivery of fertilizer N to 
farm gate accounted for 26% of the total GHG emissions, 
followed by direct and indirect emissions of 28% due to the 
application of N fertilizers to barley crop. Emissions due to 
N fertilization were 26.6 times the emission from the use of 
phosphorous, 5.2 times the emission from pesticides, and 
4.2 times the emission from various farming operations. 
Decomposition of crop residues contributed emissions of 
173 kg C0 2 eq ha -1 , or 19% of the total emission. Indian 
Head-produced barley had significantly greater grain yield, 
resulting in about 11% lower carbon footprint than Swift 
Current-produced barley (0.28 vs. 0.32 kg C0 2 eq kg -1 of 
grain). Emissions in the barley production was a linear 
function of the rate of fertilizer N applied to the previous 
oilseed crops due to increased emissions from crop residue 
decomposition coupled with higher residual soil mineral N. 
Conclusions The key to lower the carbon footprint of barley 
is to increase grain yield, make a wise choice of crop types, 
reduce N inputs to crops grown in the previous and current 
growing seasons, and improved N use efficiency. 

Keywords Carbon footprint • Life cycle assessment • N use 
efficiency • No-till practices 

1 Introduction 

Rising awareness of climate change and energy security has 
spurred the interest of quantifying carbon footprints of 
goods and services for all sectors (Wiedmann et al. 2006). 
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In agriculture, the carbon footprint of grain crops is consid¬ 
ered a new farm management indicator and can be the focal 
point in the evaluation of environmental impacts of agricul¬ 
tural activities (Weinheimer et al. 2010). Improved farming 
practices are required to minimize potential negative 
impacts of agricultural activities on environments (Liebig 
et al. 2007). In recent years, best crop management practices 
have been developed for not only producing high-quality, 
affordable food in sufficient quantity, but also reducing 
negative environmental impacts. One of the promising prac¬ 
tices is to adopt diversified cropping systems where cereals, 
oilseeds, and legume crops are arranged in well-designed 
crop sequences (Gan et al. 2010). Studies from the northern 
Great Plains of North America have shown that the diversi¬ 
fication of cropping systems can reduce production inputs 
(Zentner et al. 2001), decrease pest infestation (Krupinsky et 
al. 2002; O’Donovan et al. 1997), and lower carbon foot¬ 
prints of field crops (Gan et al. 2011c). 

Spring barley ( Hordeum vulgare L.) is a main feed grain 
crop typically grown in diversified cropping systems in 
rotation with other crops (Karamanos et al. 2010). In semi- 
arid rainfed environments, barley crop usually follows an 
oilseed crop such as flax ( Linum usitatissimum L.), sunflower 
(Helianthus annuus L.), or Brassica species such as Brassica 
juncea mustard, Brassica napus canola or Brassica rapa 
canola. The production of these oilseeds often requires high 
N fertilizer input in order to optimize seed yield and oil quality 
(Malhi et al. 2007; May et al. 2010). Under favorable growing 
conditions, canola crops respond positively to N fertilizer up 
to 180 kg Nha -1 , and some hybrid cultivars can respond to N 
fertilizer up to 250 kg Nha -1 (Brandt et al. 2007). Even in 
semiarid areas where water is the primary factor limiting crop 
production, a level of 100 kg Nha -1 of fertilization often 
results in high economic return for hybrid Brassica crops 
(Karamanos et al. 2007; Malhi et al. 2007; Cutforth et al. 
2009). 

Little has been determined whether the high level of 
fertilizer N input in the previous oilseeds would affect the 
carbon footprint of the subsequent barley crop. There exist 
huge knowledge gaps on the comparative values of carbon 
footprints of barley crops grown under various cropping 


systems. The present study may generate valuable informa¬ 
tion to begin to fill the knowledge gap. Therefore, the 
objectives of this study were to determine (1) the carbon 
footprint of barley preceded by various oilseed crops that 
were fertilized at various rates of inorganic N fertilizers the 
previous years and (2) the effect of environmental condi¬ 
tions on the carbon footprint of barley under semiarid 
rainfed conditions. We hypothesized that the type of oilseed 
crops grown the previous season and N fertilizer input to the 
previous oilseeds impacted the carbon footprint of the sub¬ 
sequent barley crop, and that the degree of the influence 
interacted by environmental conditions under which the 
crop was grown. 

2 Materials and methods 

2.1 Field experiment 

Field experiments were conducted at six environsites (i.e., 
the combinations of two locations x 3 years); they were at 
Indian Head in 2005, 2006 and 2007, and at Swift Current in 
2004, 2005 and 2006, Saskatchewan, Canada. The soil 
characteristics of the experimental locations are summarized 
in Table 1 . Treatments included the main plot, four oilseed 
species in the previous year, and the subplot, seven rates of 
N fertilizer applied to the previous oilseeds. Canola-quality 
mustard (cv. Dahinda), B. napus canola (cv. Invigor 5020), 
sunflower (cv. 63 M02), and flax (cv. CDC Bethune) were 
grown in a factorial combination with N rates of 10, 30, 70, 
90, 110, 150, and 200 kg Nha -1 , using a randomized, 
complete block design with four replicates. The oilseed 
crops were managed using recommended practices in terms 
of seeding techniques and pest control, etc. (May et al. 
2010). In the following year, spring barley (cv. Metcalf) 
was planted on the fields of standing oilseed stubbles. In 
each year, a different piece of land was used for the oilseed- 
barley sequencing plots. At Indian Head, barley was planted 
on 5 May, 6 May and 12 May, in 2005, 2006, and 2007, 
respectively, and at Swift Current on 27 April, 5 May, and 7 
May, in 2004, 2005 and 2006, respectively. A no-till drill 


Table 1 Soil characteristics of the experimental sites, Saskatchewan, Canada 


Location 

Latitude/longitude/ 

elevation 

USA soil 
description 

Canadian soil 
classification 

Texture 




Organic 

matter 

pH (water paste) 


Class 

Sand Silt 
Percent 

Clay 


Indian Head 

50°32' N; 103°40' W; 
580 m 

Typic Haplustoll 

Indian Head (heavy clay) 

Clay 

18 

39 

43 

8.3 

6.7 

Swift Current 

50° 15' N; 107°44' W; 
817 m 

Aridic Haploboroll 

Swinton (silt loam) 

Silt loam 

23 

49 

28 

3.0 

7.3 
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with hoe openers and row packers at the row spacing of 
20 cm (Indian Head) and 25 cm (Swift Current) was used for 
seeding. The plot size was 2.4 x 10.7 m at Indian Head, and 
2x14 m at Swift Current. The seeding rate was based on 
seed size, pre-seeding germination rate, and estimated field 
emergence rate of 70% to target a plant population of 120 
plants per square meter. 

At each environsite, barley crop was fertilized with urea 
which was mid-row banded 6-cm deep along with mono¬ 
ammonium phosphate (11-51-0 of N-P 2 0 5 -K 2 0), and all 
plots received the same rate of fertilizers (Table 2). The 
blanket application of fertilizers to all barley plots at each 
environsite allowed the determination of the potential effect 
of various rates of N fertilizer applied to the previous oilseed 
crops. At all sites, barley crops were grown using no-till 
management practices, where all fertilizers were applied at 
the time of sowing the crop in a single operation. Weed 
control in barley was achieved using a pre-seeding or a pre- 
emergent bum-off treatment of glyphosate at 1.2 Lha -1 of 
active ingredient (a.i.), along with bromoxynil at 500 mL ha -1 
and MCPA ester at 500 mL ha -1 plus 2,4-D 600 Amine at 
910 mL ha -1 . At physiological maturity, aboveground plant 
biomass was determined by harvesting plants in a 1.0-m area 
in each plot, and the plant samples were oven-dried at 70°C 
for 7 to 10 days and weighed. At full maturity (grain moisture 
approaching 130 gkg 1 ), entire plots were harvested with a 
plot harvester and grain dry weight determined. The N con¬ 
centration in the straw was determined using the standard 
Kjeldahl N method. Root dry weight was estimated using 
the model developed by Gan et al. (2009), where root biomass 
was a proportion of straw biomass at a given growing 
condition. 


2.2 Estimates of GHG emissions and carbon footprint 

Barley carbon footprint is largely dependent on (1) grain 
yield of the crop and (2) total greenhouse gas (GHG) emis¬ 
sions associated with the crop production, including C0 2 
emissions from energy use and N 2 0 emissions from non¬ 
energy sources. The C0 2 emission associated with urea 
fertilizers applied to the barley crop was also included. 
Using site-collected data for each of the six environsites 
coupled with empirical modeling, we estimated the GHG 
emissions deriving from (1) application of synthetic N and P 
fertilizers; (2) crop residue decomposition; (3) production, 
storage, and transportation of synthetic N and P fertilizers to 
the farm gate; (4) production of herbicides and pesticides 
and their applications; and (5) various farming operations 
such as spraying pesticides during the crop production year, 
planting and harvesting the crop. The system boundaries 
were set as the period of the life cycle from production 
inputs (such as fertilizers, pesticides, seeds), delivery of 
the inputs to farm gates, and the barley crops were harvested 
and the grain products stored in bins on the farm. Potential 
emissions along the logistics of transporting, exporting, or 
marketing the grain products leaving from the farm were not 
considered in this study because they were considered the 
same for all crop management practices evaluated in this 
study. 

It is generally known that N 2 0 as a by-product is pro¬ 
duced during nitrification and denitrification by microorgan¬ 
isms (Forster et al. 2007). The amount of direct and indirect 
N 2 0 emissions is related to the quantity of N applied along 
with environmental conditions (Dyer et al. 2010). The mag¬ 
nitude of these emissions has been studied for the Canadian 


Table 2 Monthly precipitation, 
potential evapotranspiration 

Variable 

Indian Head 



Swift Current 


(PE), fertilization, estimated 
fraction of N subject to leaching 


2005 

2006 

2007 

2004 

2005 

2006 

ana caicuiatea airect in 2 u emis- 

sion factors for synthetic N 
application and crop residue 

Precipitation May 

(mm) 

57.6 

39.0 

81.2 

83.7 

22.4 

43.5 

decomposition at Indian Head in 

June 

99.2 

80.4 

47.1 

66.2 

123.2 

99.9 

2005, 2006, and 2007 and at 

July 

59.2 

4.4 

51.4 

61.1 

21.4 

26.3 

Swift Current in 2004. 2005. and _ . . . 

2006, Saskatchewan, Canada 

August 

98.0 

11.8 

67.9 

72.3 

52.1 

24.1 

September 

4.0 

61.7 

24.4 

27.4 

40.7 

25.0 


October 

6.6 

71.7 

30.0 

21.5 

9.2 

8.8 


Total 

324.6 

269.0 

302.0 

332.2 

269.0 

227.6 


PE (mm) 

554.4 

615.9 

597.7 

520.2 

561.9 

620.7 


Fertilizer N applied to barley 
(kg Nha *) 

56 

45 

62 

35 

30 

35 


Fertilizer P applied to barley 
(kg P ha -1 ) 

25 

21 

23 

14 

17 

22 


Emission factor for N 2 0 

0.0070 

0.0050 

0.0062 

0.0092 

0.0057 

0.0033 


FRACt path 3 

0.1495 

0.1197 

0.1374 

0.1827 

0.1307 

0.0944 

a N leaching factor 
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conditions (Gregorich et al. 2005; Rochette et al. 2008). 
Using a large number of experimental observations on mea¬ 
sured N 2 0 fluxes from Canadian farmlands, Rochette et al. 
(2008) developed a simple and accurate model for deter¬ 
mining N 2 0 emission factors based on a growing season 
moisture deficit—a linear function of the ratio of growing 
season precipitation (P) to potential evapotranspiration (PE) 
as follows: 

EF = 0.022 P/ PE - 0.0048, (1) 

where EF is the emission factor with a unit of kilogram 
N 2 0-N per kilogram of N; P/PE is the ratio of P to PE 
during the growing season (1 May-30 Oct) based on long¬ 
term data (see Table 2). Soil mineral N, particularly in the 
form of nitrate, in the rooting zone is subject to leaching, 
and this N can be leached or has undergone further trans¬ 
formations to emit N 2 0. We estimated the fraction of N 
subject to leaching (FRAC LEAC h) as follows: 

FRAC LE ach = 0.3247P/PE - 0.0247 (2) 

Using the method developed by the Intergovernmental 
Panel on Climate Change (IPCC) adopted for Canadian 
conditions (IPCC 2006), emissions of N 2 0 from synthetic 
N applications were estimated as follows: 

CO 2 eq SNF _ N20 = gsNF x {(FRACqasm x EF vd ) + EF 

44 

+(FRACleach x EFleach)} x x ^98 

( 3 ) 

where C0 2 eq S NF-N2o was the total emissions from the syn¬ 
thetic N fertilizer application (kilogram C0 2 eq per hectare), 
6snf was the quantity of synthetic N fertilizer applied 
(kilogram N per hectare), FRAC gasm was the fraction of 
synthetic N fertilizer that volatilized as NH 3 ~ and NO x -N 
(FRAC gasm = 0.1 kg Nkg -1 N), EF VD was the N 2 0 emis¬ 
sion factor for volatilized NH 3 ~ and NO x -N (EF VD = 
0.01 kg Nkg -1 N), EFeeach was the N 2 0 emission factor 
for nitrate leaching (EF L each= 0.0075 kg Nkg -1 N), 44/28 
was the conversion coefficient from N 2 0-N to N 2 0, and 298 
was the global warming potential of N 2 0 over 100 years 
(IPCC 2006). In the case where urea was applied to the 
barley crop, the carbon contained in the urea is often 
released as C0 2 during hydrolysis (IPCC 2006). Thus, 
emissions of C0 2 from urea-based N fertilizer were calculated 
as follows: 

12 44 

C0 2 eq SNF _ C02 — gsNF-uREA x 28 X L2 ’ 4 

where C0 2 eq SNF .co2 was the emissions of C0 2 from the urea 
application (kilogram C0 2 eq per hectare), 2sNF-urea was the 
quantity of urea-based N fertilizer applied (kilogram N per 
hectare), 12/28 was the ratio of C to N, and 44/12 was 


converted C to C0 2 . When the barley crop was harvested, 
straw and roots were returned back to the soil. The N 
contained in the barley crop residues provided additional 
source of N for nitrification and denitrification, leading to 
N 2 0 emissions. The quantity of barley crop residue N 
(Gcrd) was obtained using the aboveground and below¬ 
ground crop residue biomass multiplied by its respective N 
concentration. Thus, emissions from barley crop residue 
decomposition were calculated as follows: 

C0 2 eq CRD = 2 crd x {EF + (FRACleach x EFleach)} 

44 

*28 * 298 ( 5 ) 

The Haber-Bosch process that converts N 2 and H 2 gases 
into NH 3 is known to be energy and emission intensive (Gan 
et al. 2011a). Based on available data from the studies 
conducted in North America (Lai 2004; Rochette et al. 
2008), we estimated emissions from the production, trans¬ 
portation, storage, and delivery of fertilizers to farm fields 
using an emission factor of 4.8 kg C0 2 eq kg -1 of N and 
0.73 kg C0 2 eq kg -1 of P 2 0 5 , multiplied with the amount of 
N and P fertilizers on a per-hectare basis. 

Herbicides and fungicides were used in the produc¬ 
tion of barley crops. Although emission factors for each 
individual pesticide are not available at the present time, 
we assumed that the emission during the process of 
production, transportation, storage, and field application 
was similar among pesticides within a similar category. 
Thus, emissions due to pesticides were estimated based 
on active ingredient of a fungicide or herbicide product, 
with the emission factor of 23.1 kg C0 2 eq kg -1 of a.i. 
for herbicides and 14.3 kg C0 2 eq kg -1 of a.i. for fungicides. 
The emissions associated with various farming operations 
such as sowing, spraying, and harvesting were estimated using 
a factor of 14 kg C0 2 eq ha -1 for no-till planting, 
5 kg C0 2 eq ha -1 for herbicide and fungicide spraying, and 
37 kg C0 2 eq ha -1 for harvesting (Lai 2004; Gan et al. 2011b). 

With the above factors, we estimated (1) total GHG 
emissions per unit of areas for the production of barley, 
expressed as kilogram C0 2 eq per hectare and (2) carbon 
footprint as emissions per kilogram of barley grain produced 
under the specific growing conditions, expressed as kilo¬ 
gram C0 2 eq per kilogram of grain. 

2.3 Statistical analysis 

Data were analyzed using the PROC MIXED model of SAS 
(Littell et al. 1996), where oilseed crop type of the previous 
year and fertilizer N rate applied to the previous oilseeds 
were designated as fixed effects and replicates as random 
effects. In the analysis, N fertilizer was considered a class 
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variable, whereas the various N rates received by previous 
oilseeds was considered a continuous variable. Therefore, 
all interactive responses of crop types and environsites to the 
various N rates were determined by performing linear and 
nonlinear regressions in the analysis of variance. Following 
the ANOVA analysis, covariance analysis was performed to 
determine the possible main factors contributing to the var¬ 
iance observed in emissions and carbon footprint of barley. 
Residual soil mineral N in the 0-60-cm depth prior to 
seeding was treated as a covariable in the analysis. In all 
analyses, significant effects were declared at ^<0.05. 


3 Results 

3.1 Productivity of previous oilseeds and re-cropped barley 

Detailed results on the biomass, seed yield, and quality of 
the four oilseed species grown the previous year have been 
discussed in a previous publication (May et al. 2010). To 
facilitate the carbon footprint estimate of the subsequent 
barley crop, we presented a summary of the productivity 
of oilseed crops. In brief, all four oilseed species had a 
curvilinear increase in seed yield as fertilizer N rates increased 
from 0 to 200 kg Nha -1 , with the largest yield response 
observed in B. napus canola to the high end of N rates. The 
majority of the increase in flax seed yield occurred as the N 
rate increased from 10 to 90 kg ha -1 , while most of the 
increase in seed yield of B. juncea canola and sunflower 
occurred as N rates increased from 10 to 70 kg ha -1 . For all 
oilseeds tested, their protein concentration increased but oil 
concentration decreased as the N rate increased. With the 
current available cultivars in each species, B. napus hybrid 
canola produced the greatest oilseed yield and required the 
highest N input. 

The grain yield of barley crops differed significantly 
among environsites (Fig. la). Barley grown at Indian Head 
in 2005 produced the highest yield (4,682 kg ha -1 ), whereas 
the barley at Swift Current in 2005 yielded the lowest 
(1,368 kg ha^ 1 ). Barley grain yield was also influenced 
significantly by crop species grown the previous year (see 
Fig. la), but the majority of the differences was due to 
sunflower crops which resulted in the lowest grain yield of 
barley in most environsites, except at Swift Current in 2006 
where no differences were found. Overall, the barley pre¬ 
ceded by flax, canola, and juncea crops produced similar 
grain yields at four of six environsites, with the barley 
preceded by canola at Indian Head in 2006 and Swift 
Current in 2005 having greater grain yields than being 
preceded by the other oilseed species. 

The grain yield of barley crops was linearly related to the 
amount of fertilizer N applied to the previous oilseeds (see 
Fig. lb). Increasing rate of fertilizer N applied to the 



U) 



Canola 

Flax 

Juncea 

Sunflower 








IH-2007 SC-2004 SC-2005 SC-2006 


2 Environsites 

0 ) 



Fertilizer N (kg N ha' 1 ) 
applied to the previous oilseed crops 


Fig. 1 Average grain yield of barley crop a that was preceded by 
different oilseed crops (across all N rates) at Indian Head (. IH) in 2005, 
2006, and 2007 and at Swift Current (SC) in 2004, 2005, and 2006 and 
b in relation to the rate of N fertilizer applied to the previous four 
oilseed crops, in Saskatchewan, Canada. Bars are standard errors 


previous oilseed crops significantly increased the grain yield 
of succeeding barley crops at all environsites except at Swiff 
Current in 2005. The magnitude of the yield response was 
interacted with environsites; every 1 kg of fertilizer N applied 
to the previous oilseeds increased barley grain yield by 5.6 kg 
at Indian Head in 2005,4.4 kg at Swift Current in 2004, 3.9 kg 
at Indian Head in 2007, 2.0 kg at Swift Current in 2006, and 
1.4 kg at Indian Head in 2006, but decreased by 1.5 kg at Swift 
Current in 2005. 

3.2 Main contributors of emission estimates in barley 

The long-term (1971-2010) average growing season (1 
May-31 Oct) precipitation was 264 mm at Swift Current 
and 317 mm at Indian Head, and PE was 635 mm at Swift 
Current and 605 mm at Indian Head. Using Eq. 1, we 
estimated that direct emission factors for soil N 2 0 from 
synthetic N application and crop residue decomposition 
were 0.0067 kg N 2 0-N kg -1 N at Indian Head and 
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0.0044 kg N 2 0-N kg -1 N at Swift Current. On a relative 
basis, the direct emission factor at Indian Head was about 
50% greater than that at Swift Current. However, weather 
conditions varied widely among environsites during the 
course of the experiment (see Table 2). The most noticeable 
is that at Indian Head, the growing season precipitation was 
higher, but the PE was lower than the long-term average in 
2005. The opposite was true in 2006 while the 2007 grow¬ 
ing season precipitation and PE were near the long-term 
average. At Swift Current, the growing season precipitation 
was greater, but the PE was lower than the long-term 
average in 2004. As a result, the estimated direct N 2 0 
emission factor for 2004 at Swift Current was 0.0092 kg 
N 2 0-N kg -1 N, the highest among all environsites. In 2006, 
lower than the long-term average growing season precipita¬ 
tion and average PE at Swift Current resulted in the lowest 
direct N 2 0 emission factor (0.0033 kg N 2 0-N kg -1 N) among 
all environsites. 

Many factors contributed to the GHG emissions in the 
production of barley crops. Covariance analysis revealed 
that 89% of variance in GHG emissions was due to growing 
environments (Table 3). Nitrogen fertilizer applied to the 
previous oilseeds contributed 7% of the variation in the total 
GHG emissions, and the other factors, albeit being impor¬ 
tant to barley production, contributed little to the variation 
of total emissions. Overall, environsites had a highly signif¬ 
icant effect on total emission (see Table 3). Averaged across 
the four oilseed species grown the previous year, the suc¬ 
ceeding barley crop grown at Indian Head emitted GHGs of 
1,003 kg C0 2 eq ha -1 , or 53% greater than the emission of 
barley crop grown at the drier Swift Current site (Table 4). 
Even at the same experimental site, total emissions differed 
significantly between years mainly due to precipitation and 


evapotranspiration. At Indian Head, the 2005 barley crop 
had total emissions 52% greater than the 2006 barley crop, 
whereas at Swift Current, the 2006 barley crop had total 
emission 32% less than the 2004 barley crop. 

Nitrogen, phosphorous, and pesticides are the main 
inputs in the production of barley crops. Estimated 
emissions from the production, transportation, storage, 
and delivery of inorganic fertilizers to the farm gate 
accounted for an average of 26% of the total emissions 
(see Table 4), and the direct and indirect emissions from 
the application of N fertilizer to the barley crop added 
another 28%. Thus, the production and application of 
synthetic N together contributed 53% of the total emis¬ 
sions. Such a proportion (i.e., fertilizer-N-related emis¬ 
sions in the total emission) was nearly the same across 
all environsites, even though the absolute emission val¬ 
ues varied largely among environsites. For example, the 
emission due to the application of N fertilizer was 
129 kg C0 2 eq ha -1 at Swift Current in 2006 but was 
350 kg C0 2 eq ha -1 at Indian Head in 2007. Averaged 
across the six environsites, the GHG emissions due to 
the use of fertilizer-N was 26.6 times the emission from 
the use of phosphorous, 5.2 times the emission from 
pesticide supply, and 4.2 times the emission from vari¬ 
ous farming operations. 

The decomposition of barley crop residues contributed 
direct and indirect emissions, averaging 173 kg C0 2 eq ha -1 , 
or 19% of the total emission (see Table 4). This portion of the 
emission varied among environsites, largely due to the varia¬ 
tion in the amount of barley straw and root biomass and the 
difference in P/PE ratio. For example, the barley grown at 
Indian Head in 2005 produced 6,720 kg ha -1 of straw and 
2,620 kg ha -1 of roots, resulting in the highest emission from 


Table 3 Summary of covari¬ 
ance of analysis for total emis¬ 
sions and carbon footprint for 
barley grown at various previous 
rates of N fertilizer at eight 
environmental sites (location- 
year combinations) in Saskatch¬ 
ewan, Canada 


Bold numbers represented that 
the effect was significant at 
P<0.05 

Tn the covariance analysis, N 
rate was considered as a contin¬ 
uous variable, whereas nitrogen 
was consider as a ‘class’ factor 

b Env environsites, C preceding 
crop type, N N rates in previous 
crops 


Source 

DF 

Total GHG emission in the 
production of barley crops 

Carbon footprint of barley 
crops 

Sum of squares 

% variance 

Sum of squares 

% variance 

Environsites 

5 

8,656,697 

88.7 

1.0750 

72.0 

Preceding crop type 

18 

124,322 

1.3 

0.1741 

11.7 

N rates in previous crops a 

1 

669,688 

6.9 

0.0324 

2.2 

N x N 

1 

38,028 

0.4 

0.0068 

0.5 

Nitrogen a 

5 

1,678 

0.0 

0.0046 

0.3 

C x Env b 

5 

95,472 

1.0 

0.1237 

8.3 

N x N x Env 

5 

13,259 

0.1 

0.0095 

0.6 

Env x nitrogen 

25 

9,192 

0.1 

0.0068 

0.5 

Seed yield of barley 

1 

57,552 

0.6 

0.0123 

0.8 

Residual soil mineral N at seeding 

1 

89,863 

0.9 

0.0195 

1.3 

Crop residues of barley 

1 

1,240 

0.0 

0.0019 

0.1 

Residual 

122 

0 

0.0 

0.0264 

1.8 


4^ Springer 








Int J Life Cycle Assess (2012) 17:635-645 


641 


23 

£ 

o 
+-> 

<D 

CO 

C 

o 

C/3 

c/3 

2 

CD 

C 

c3 

co 

C 

O 

c/3 

c/3 

2 

1) 

^2 


c3 

c3 

Ch 

c3 

u 


xT § 
13 ^ 

• 1 —I QD 

>. X 

.S B 

s 

bD g 
g C 

Vh • 

CD 

3 


CD 

3 

f3 


£ 

2 

60 

co 

a 

o 

o 


cr 

4 _ U 

7h ^ 

?h .2 O xj 

O b r ) ,— s OD 

x ^ cd 

c3 o to I ft “ 

r \ o Xs3 «0 ^-H 

U « o 



C/3 

C /3 

• rH 
2 
w 


o 

.*H 

<H-H 

£ 

o 

• rH 

C /3 

C /3 


§ 


I'ss 

3 >,3 

b b? o 

H T3 
g O Ph OO 
r « h d d 

W a m a 


2 

c§ 

<H-H 

a 

o 

• rH 

C /3 

C /3 


<D 

•§ 


Ph 

o 


to 
O 

2 
a p 


co 

CD 

pH 


2 g cd 
w o ^ 


o 

CD 


2 

eg 

G—I 

a 

o 


a 

o 


■ „■ o 

co i-j 
co H 


a 

o 

• rH 
H-H 

ctf 

o 

* rH 

'Hh 

& 


• h -d " 7 

2 o ^ £ 

r « p P 'w 
W pH Cb O 


Pt 

O 

• t-H 
4-> 

cb 

CD 


o 

c£h Ph 
, CD 
PP N 

o x 

• .-H . t-H 

co -r 
•r X H-H 

S b g 

W £h cb 


o 

^ ^ g 

pp N .2 

o b -P 

•r rpp o 

m -h b 

CO r X 

CD X 
^ O 

w £ a 


pp 

o 

• rH 

1/2 r^H 

co CT 
■pp tD 

CD O 

13 O - 

-K 60 I 
£ <* .§ 


I^>rH 

23 1 

P ^ b rP 

§ 3 £ gp 
O o ^ 


cb 

£ 


PP 

o 

• t-H 
H—' 

a 

CD 

o 


CP 

CD 

CD 

Ph 

CD 

Ph 


f cd 

M-H 

i/i 

a 




C/3 

o 

• ^H 

H—' 


cr 


r—* 

a 


CD 


• rH 
$-H 

$H 

<D 

On 

bti 

CN 

o 

1 

c3 

> 

O 

^4 

o 



PP 

CD 

CD 

Ph 

CD 

Ph 


cr 

CD 

CN ’—i 

60 Q 1 c3 
X O X 


cP 

CD 

CD 

pH 

CD 

Ph 


cr 

CD 

(N^h 

bO Q 1 cb 
X U X 


PP 

CD 

CD 

Ph 

CD 

Ph 


SP 

CD 

CD 

Ph 

CD 

Ph 


X 

CD 

CN 

o 

u 


60 'cb 

X X 


pp 

CD 

CD 

Ph 

CD 

Ph 


ct 

CD 

<N^H 

bO Q 1 cb 
X U X 


cn ^h in 
cn oc — 
cn cn co CN 

o o o o 


o\ 


ON 


cn 

r- 


cn 

o 

o 

CN 

nd 

cd 

CD 

ffi 

§ 
• rH 

r rd 

pp 


X CN cn 


— 1 CN CN — 
O CN CN O 


T vf O (N 


O IN OO CN 

-H o o o 


ON X G> IN 
CN t— 11 i—i CN 


N i—i X 

^Ct- CN O CN 
cn — 1 CN CN 


cn nf- On On 
cn, m cn n- 


O X 
cn cn 


cn 


ON 

CN 


o cn cn o 
cn O on nt - 
cn cn cn cn 


cn oo oo o 
CN CN CN CN 


IN in NO NO 

NO ih ON no 
CN CN CN '— 1 


ON 

H- ^ 

N O 

I ' l-H 


m- 

cn 

oo 


cd 

cn cn -^d- 
oo n cn n 

x N N" nj 

cn cn rs rs 

nr cn cn cn 


X N N" 

o o o 
o o o 

CN CN CN 


X oo 
CN NO 

m- cn 
o o 


IN IN 


nf- 

ON ON 


OO ON 


'vt oo 
o o 


m cn 


CN IN 
IN X 


CN CN 
in in 


IN cn 
CN CN 


O ON 

un cn 


un 

cn 


m 

m- 


o 

cn 


IN 

NO 


cn un 
x x 
un un 


OO ON 
NO *“H 
cn *—i 

C\ CN 

^H <N 

m x 
o o 
o o 

CN CN 


PP 

CD 

h_ Ip 

3 B 

* u 
GO 


170 
H—> 

PP 

0D 


cb 

CD 

S—i 


cb 

170 

170 

O 

S—i 

o 

cb 

nd 

OD 

bD 

cb 

Ph 

OD 


the crop residue decomposition at 347 kg C0 2 eq ha -1 , or 5.1 
times the same category of emissions by the barley crop 
grown at Swift Current in 2006. 

Emissions associated with pesticide use and various 
farming operations in the production of barley crops were 
lower than the emission from the use of N fertilizers (see 
Table 4). The barley crop received a pre-seeding bum-off 
treatment with glyphosate plus two in-crop applications of 
herbicides, one for grassy weed control and the other for 
broadleaf weed control. 

3.3 Carbon footprint versus preceding crops and N rates 

Covariance analysis revealed that the carbon footprint of 
barley crop was significantly influenced by the choice of 
oilseeds grown in the previous year (see Table 3). At Swift 
Current in 2005, when the barley was preceded by sunflower, 
it had a significantly higher carbon footprint than when pre¬ 
ceded by other oilseed crops (Fig. 2). Sunflower is a deep¬ 
rooting crop which depleted soil water throughout the profile. 
As a result, barley grown after sunflower decreased grain yield 
by an average of 23% compared with barley grown after flax, 
37% decrease compared with barley after B. juncea canola, 
and 59% decrease compared with barley after B. napus canola 
(Fig. la). The significant effect of previous oilseed types on 
the grain yield of re-cropped barley led to the significant 
differences in its carbon footprint. 

Both emissions and carbon footprint of barley crop also 
were significantly influenced by the rate of N fertilizer 
applied to the previous oilseeds (see Table 3). More N 
fertilizer applied to the previous oilseeds gave rise to greater 
GHG emissions in the subsequent barley crop. At all envi- 
ronsites, the amount of residual soil mineral N measured 
prior to seeding the barley crop was increased as the amount 
of N applied to the previous oilseeds more than 90 kg Nha -1 



Environsites 

Fig. 2 Carbon footprint of barley crop preceded by different oilseed 
crops (averaged across N rates) at six different environmental sites in 
Saskatchewan, Canada. Bars are standard errors 
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(Fig. 3a). This effect was consistent for all the oilseed 
species evaluated (see Fig. 3b). Regression analysis showed 
that total emission in the production of barley crop was a 
linear function of the rate of fertilizer N applied to the 
previous oilseeds at four of the six environsites (Fig. 4a). 
At those environsites, increased N fertilizer application to 
the previous oilseeds increased the total emission of the 
subsequent barley crops. This was largely due to the greater 
amounts of straw and roots (data not shown) produced by 
barley crops grown in the high-N soil, causing more GHG 
emission from the decomposition of the barley resides. The 
relationship between N rates applied to the previous oilseed 
crops and the emission of the subsequent barley crops was 
consistent across various environsites (see Fig. 4a). 

3.4 Carbon footprint versus environments 

Covariance analysis revealed that the carbon footprint of 
barley differed significantly among environsites; 72% of 
the variation in the barley carbon footprint was due to 
environments (see Table 3). The barley grain produced at 
Swiff Current in 2005 had the highest carbon footprint at 
0.426 kg C0 2 eq kg^ 1 of grain, 70% greater than the carbon 
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Fig. 3 The relationship between the rates of fertilizer N applied to the 
previous oilseed crops and the residual soil mineral N (NH 4 + and 
N0 3 ~) in the 0-60-cm depth measured prior to the seeding of barley 
a at various site-years and b on the fields with different oilseed crops 
the previous year 
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Fig. 4 The effect of the rate of fertilizer N applied to the previous 
oilseed crops on a the greenhouse gas emissions and b the carbon 
footprint of the subsequent barley crop at Indian Head ( IH) in 2005- 
2007 and at Swift Current ( SC) in 2004-2006. Bars are standard errors 


footprint value of the barley produced at Indian head in 
2005 (0.252 kg C0 2 eq kg -1 of grain) (see Table 4). On 
average, Swift Current-produced barley had about 11% 
greater carbon footprint than Indian Head-produced barley 
(0.317 vs. 0.281 kg C0 2 eq kg -1 of grain). Although total 
emissions associated with the barley production at Indian 
Head were generally greater than the emission at Swift 
Current, the barley produced at Indian Head had significant¬ 
ly greater grain yield. Carbon footprint was a function of 
grain yield and total GHG emissions. Increasing grain yield 
decreased carbon footprint in barley, and the magnitude of 
the decrease depended on total GHG emissions. 

The relationship between N rate received by the previous 
oilseed crops and the carbon footprint of the subsequent 
barley crops was significantly altered by environments (see 
Fig. 4b). At Swiff Current in 2005, the increased N rates to 
the previous oilseeds significantly increased barley carbon 
footprint; the opposite was true at Indian Head in 2006 when 
the increased N rates to the previous oilseeds significantly 
decreased barley carbon footprint; at the other environsites, 
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the barley carbon footprint did not change with the N rate of 
the previous oilseeds. The large variation in barley carbon 
footprint between environsites was a result of the different 
magnitudes of responses of grain yield (see Fig. lb) and the 
total emission (see Fig. 4a) to the rate of N applied to the 
previous oilseed crops. At a given N rate, the barley grown 
at Swift Current in 2005 had its carbon footprint significantly 
greater than the values at the other environsites. In that partic¬ 
ular year, high rainfall in the month of June (see Table 2) 
caused excessive vegetative growth in barley, but severe 
drought occurred in July prevented the photo synthetic materi¬ 
als from being remobilized into grain sinks, resulting in low 
seed yield (see Fig. 1). 

4 Discussion 

4.1 Estimation methodology and main contributors to barley 
carbon footprint 

Total GHG emissions associated with the production of 
barley crops included C0 2 emissions from energy use and 
N 2 0 emissions from non-energy sources. To estimate the 
total emission accurately, we measured various emission 
contributors from the site-specific plot experiment with an 
approach similar to that used previously (Rochette et al. 
2008; Gan et al. 2011b). Also, we employed the IPCC 
methodology (IPCC 2006) adopted for Canadian conditions 
(Rochette et al. 2008) to establish the emission factors of 
each emission contributor. We feel that the approach of 
taking the site-specific measurement combined with adopt¬ 
ing empirical relationships is quite effective in evaluating 
the impact of farming practices on the carbon footprints of 
field crops. 

Synthetic N fertilizers used in the production of barley 
crops contributed the greatest percentage (53%) of the total 
emissions. Of which about 28% of the emission came from 
direct and indirect emissions of N 2 0 from the application of 
N fertilizers in the field, and the other 26% coming from the 
production, transportation, storage, and delivery of N fertil¬ 
izers to farm gates prior to farm use. When barley crop was 
harvested, the straw and roots were returned back to the 
field under no-till management practices. These plant litters 
served as an important N source for nitrification and deni¬ 
trification, contributing directly and indirectly to N 2 0 emis¬ 
sions. In our study, the decomposition of barley straw and 
roots contributed an average of 19% to the total emission. In 
a previous study, Gan et al. (2011c) reported that straw and 
roots in durum wheat contributed about 25% of its carbon 
footprint, whereas the decomposition of oilseed straw and 
roots contributed about 10% of the total carbon footprint 
(Gan et al. 2011b). The large difference between oilseeds 
and cereals is because the cereals have greater amounts of 


straw and root C (typically >2,200 kg Cha -1 ) than oilseeds 
(typically <1,700 kg Cha -1 ) under similar growing condi¬ 
tions. Also, the magnitude of the contribution from the 
decomposition of crop residues is related to the net produc¬ 
tivity of a crop (Forster et al. 2007), N concentrations of the 
plant residue (Janzen et al. 2006), and growing conditions 
such as soil moisture and temperature (Flynn et al. 2005; 
Merrill et al. 2007). 

4.2 Barley carbon footprint versus cropping practices 

Under a specific environment, the GHG emission in the 
production of barley crop was influenced by preceding crop 
species and N fertilization applied to the previous crops. 
Sunflower grown in the semiarid area depleted soil moisture 
in the deeper layers (Angadi and Entz 2002), decreased 
grain yield of the subsequent barley crops, and increased 
barley carbon footprint. At most environsites evaluated, the 
total emission of barley crops was linearly associated with 
the amount of N fertilizer applied to the previous oilseeds. 
The higher the N fertilizer applied to the previous crops, the 
greater the emission the subsequent barley crop. Fertilizer N 
applied to the previous oilseeds at the rate greater than 
90 kg Nha -1 significantly increased the GHG emission of 
the succeeding barley. The oilseed crops grown the previous 
year had little or no response to fertilizer N greater than 
90 kg Nha -1 (May et al. 2010), and the left, unused N 
increased the amount of residual soil mineral N significantly 
(see Fig. 3). The residual soil mineral N had a positive effect 
on barley straw and root biomass and thus increasing total 
emissions. 

However, higher GHG emissions in barley production 
due to N fertilization or preceding crops did not always 
mean a higher carbon footprint, depending on whether or 
not the cropping practices resulted in greater grain yields. 
Increased grain yield decreased the carbon footprint of bar¬ 
ley significantly, and the magnitude of the change in carbon 
footprint depended on total GHG emissions. This was the 
case in our study that at Swift Current in 2005, the increased 
N fertilization to the previous oilseeds significantly in¬ 
creased GHG emission, but the grain yield of the suc¬ 
ceeding barley crop was negatively affected by drought 
which occurred during the reproductive growth period, result¬ 
ing in the high carbon footprint of barley. In contrast, at Indian 
Head in 2006, the increased N fertilization to the previous 
oilseeds significantly increased both GHG emission and grain 
yield, but the yield increase was more substantial than the 
GHG increase, thus lowered the carbon footprint of barley. 
Our finding was in agreement with van Groenigen et al. (201 0) 
who used the term yield-scaled N 2 0 emissions as N 2 0 emis¬ 
sions per unit of grain yield. Our study indicates that the key to 
lower the carbon footprint of barley crops is to increase grain 
yield and improve N use efficiency. It can be reasonably 
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expected that all cropping practices that help increase grain 
yields and improve N use efficiency will lower the carbon 
footprint of the crop. 

4.3 Potential limitations 

In barley production, the most significant emission related to 
the energy use is the production of synthetic N fertilizers, 
accounting for about 28% of the total emissions. The produc¬ 
tion of pesticides contributed a similar amount of emissions as 
various farming operations, each accounting for about 14% of 
the total emission. Thus, more than 53% of the total emissions 
for barley production were directly related to the on-farm and 
off-farm energy use. In contrast, emissions due to fossil C0 2 
was very low in barley production, and the potential impacts 
of cropping practices on fossil C0 2 emissions were assumed 
the same in the present study because the treatments evaluated 
in our plot study did not cause large differences in fossil C0 2 
emissions. However, there are other terms in the farm energy 
budget such as gasoline consumption for on-farm/farm-owned 
transport, manufacture and supply of farm machinery, heating 
fuels (natural gas), and farm electrical power, among others. 
We were unable to measure those energy terms due to the 
small scale of the plot study, therefore, the estimated carbon 
footprint values mainly showed how cropping management 
practices affected the carbon footprint of barley and did not 
intend to be a complete value. Also, the exact amount of 
emissions might change slightly if those un-measured energy 
terms were considered, although adding those on-farm fuel 
uses will unlikely change the overall conclusions of the pres¬ 
ent study. 

The method used in estimating emissions of N 2 0 from 
synthetic N application and crop residue decomposition in 
the study was primarily based on the methodology of 
Rochette et al. (2008) who used many years of field flux 
measurements of N 2 0 at Swift Current, the site where the 
present study was conducted. In our calculation, a generic 
NH 3 release rate was used according to IPCC recommenda¬ 
tion although the NH 3 emissions may differ between N 
fertilizers, and indirect emissions (such as N 2 0 from NH 3 
or C0 2 from CO) were not taken into account. Also, carbon 
change in the soil due to crop residues was not considered in 
our calculation of carbon footprint; this followed the guide¬ 
line of PAS 2050 (BSI and Carbon Trust 2011) where it is 
indicated that the assessment of emission must include GHG 
emissions arising from direct land use change but must not 
include those arising from soil carbon changes in existing 
agricultural systems. 

Furthermore, the relationships among soil and climatic 
environments, as agronomic practices and N 2 0 emissions 
are complex (Osborne et al. 2010). Emissions of C0 2 , N 2 0, 
and CH 4 are often complicated by crop, water management, 
biogeochemical characteristics of soil (i.e., pH, salinity, C/N 


ratio, bulk density, redox potential), microbial community, 
weather, soil fertility status, and crop management practices. 
These diverse factors and their relationships to gas produc¬ 
tions make it complex and difficult to generalize. While the 
model and the IPCC method to estimate C0 2 eq were tested 
on many Canadian farms, these models might not the best 
describing seasonal and temporal variations of GHG emis¬ 
sions in barley production. Nevertheless, we believe that the 
findings on the effect of the specific cropping practices (the 
choice of oilseed species and the rate of N fertilization 
applied to the oilseeds) on the carbon footprint of barley 
was novel and valuable for land management and modelers 
who are interested in lowering the carbon footprints of field- 
grown crops. 

5 Conclusions 

More than 70% of the variation in the estimated carbon 
footprint of barley was caused by environments. On average, 
Indian Head-produced barley had a carbon footprint of 
0.281 kg C0 2 eq kg -1 of grain, significantly lower (11%) than 
Swift Current-produced barley (0.317 kg C0 2 eq kg -1 of 
grain); this was largely because the barley produced at 
Indian Head had significantly higher grain yield than the same 
crop produced at Swiff Current. Nitrogen fertilization in the 
barley production contributed more than 50% of the total 
GHG emission. Increasing N fertilizer to the previous oilseed 
crops significantly increased the GHG emission of the suc¬ 
ceeding barley crops, suggesting that it is possible to shift 
emissions from one crop to another within a crop rotation 
sequence. The carbon footprint of barley was a function of 
grain yield and total GHG emissions; increased grain yield 
decreased the carbon footprint, and the magnitude of the 
change depended on the value of total GHG emissions. The 
key to lower the carbon footprint of barley is to increase barley 
grain yield through best crop management practices such as 
optimizing N fertilizer input and maximizing N use efficiency. 
The present study had the system boundary set from the 
production of inputs, the delivery of inputs to farm gates, 
and the barley grain harvested. A full life cycle assessment 
may be required to estimate carbon footprints for a whole 
production-marketing chain or under various choices of 
boundaries; this can be achievable in the near future as more 
advanced scientific tools are becoming available. 
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